Spin current and polarization in impure 2D electron systems with spin-orbit coupling 
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We derive the transport equations for two-dimensional electron systems with spin-orbit interaction 
and short-range spin-independent disorder. In the limit of slow spatial variations of the electron 
distribution we obtain coupled diffusion equations for the electron density and spin. Using these 
equations we calculate electric-field induced spin accumulation and spin current in a finite-size 
sample for arbitrary ratio between spin-orbit energy splitting A and elastic scattering rate r~^. We 
demonstrate that the spin-Hall conductivity vanishes in an infinite system independent of this ratio. 
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Introduction. The subject of the novel and quickly 
developing field of spintronics is the transport of elec- 
tronic spins in low-dimensional and nanoscale systems. 
A possibility of coherent spin manipulation represents 
an ultimate goal of this field. Typically, spin transport is 
strongly affected by a coupling of spin and orbital degrees 
of freedom. The influence of the spin-orbit interaction 
is two-fold. The momentum relaxation due to diffusive 
scattering of carriers, e.g. by disorder, inevitably leads 
to a spin relaxation and destroys spin coherence. On 
the other hand, the controlled orbital motion of carri- 
ers can result in a coherent motion of their spins. Thus, 
spin-orbit coupling is envisaged as a possible tool for spin 
control in electronic devices. In particular, it is possible 
to generate spin polarization and spin currents by apply- 
ing electric field, the phenomenon known as the spin-Hall 
effect. 

Although the study of spin-Hall effect recently evolved 
into a subject of intense research [ll l^j Id, ljj a IS. tl 
y, y, llll lul vM , the issue remains highly controversial. 
Sinova et al. l2!| have predicted that in a clean, infinite, 
homogeneous 2DES the spin-current j^ = \{(Ti,Vk} de- 
velops a non-zero expectation value under an external 
electric field E. (Here ^<t and v are the operators of the 
electron spin and velocity, respectively.) The spin-Hall 
conductivity, defined as the ratio Gsh = ^Jy/E^, was 
predicted to have a universal value UsH — ^, indepen- 
dent of the magnitude of the spin-orbit energy splitting 
A. The effect of impurity scattering on a spin current 
has been discussed in Refs. |j,J2LUii lUl- Schhemann 
and Loss |j| and Burkov et al. [llj found that the spin- 
Hall conductivity disappears in the dirty limit A ^ t^^, 
reaching the universal value only for sufficiently clean 
regime, A 3> t~^. The clean regime has been analyzed 
by Inoue et al. |9j , who argued that the spin-current com- 
pletely disappears due to vertex corrections. Recently, 
Dimitrova 12J obtained the universal value independent 
of the relation between the spin-orbit splitting A and the 
impurity scattering rate. 

Because the spin-current is not measurable directly, its 
physical meaning is obscure. In the presence of spin-orbit 
interaction, electron spin is not a conserved quantity, and 



a spin current is not directly related to the transport 
of spins. In particular, Rashba [3 demonstrated that 
spin current can be non-zero even in equilibrium, as the 
symmetry of an isotropic spin-orbit Hamiltonian allows 
non-zero in-plane currents jy — —j}j, 7^ 0. A more mean- 
ingful quantity is spin polarization (spin accumulation) 
rather than a spin current. Equilibrium currents do not 
lead to spin-accumulation. It remains unclear whether 
the predicted nonequilihrium spin-Hall currents j^ accu- 
mulate near sample boundaries. Bulk polarization has 
been studied in both the three-dimensional [13 and two- 
dimensional [iJI electron systems in the electric field. 

In this Letter, we develop a consistent microscopic ap- 
proach to spin transport in impure 2DES. We derive a 
quantum kinetic equation which describes the evolution 
of a density matrix of a non-interacting 2DES. For length 
scales exceeding the mean free path, this equation re- 
duces to a modified diffusion equation. We then compute 
spin polarization and spin current in a general situation 
when the finite-size system is driven out of equilibrium 
by an external electric field as well as by the density gra- 
dient. We find that the spin current actually vanishes in 
an infinite system for arbitrary Ar. 

However, in a mesoscopic conductor connected to two 
massive metallic contacts, non-equilibrium spin currents 
jy flow in the vicinity of the contacts (as shown in Fig. 1). 
A non-zero spin-Hall effect can also be achieved in an in- 
finite system by applying a finite frequency electric field. 
We evaluate the ac spin-Hall conductivity, which is max- 
imal for a frequency of order of the spin relaxation rate. 
This result is instructive for making a connection with 
previous works and clarifying the 'universality' issue of 
the spin-Hall conductivity. 

Kinetic equation. Non-interacting electrons in an 
asymmetric quantum well can be described by a single 
particle Hamiltonian 
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where p = —iK^ is electron momentum, m is the effec- 
tive mass, A(t) is a vector potential of the uniform elec- 
tric field E — —A/c, and 17 is proportional to the electron 



spin operator. (We neglect terms qubic in p.) The dis- 
order potential Ui is assumed to be random, short-range, 
and spin- independent. For the isotropic ( "Rashba" ) spin- 
orbit interaction |l5|, fj = z x a, where o" are the Pauli 
matrices. To describe a non-equilibrium state of the sys- 
tem, we use the Keldysh approach [ig . We introduce the 
retarded and advanced Green's functions G^ and G^, 
and Keldysh function G^ satisfying Dyson's equation 
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Here the lower bar denotes a matrix in Keldysh space. 
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is a density of states per 



spin direction. Neglecting weak-localization effects, one 
can relate the self-energy S to the Green's function G by 
a standard disorder averaging technique |l7| . 
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We consider only the limit where r ^ and A are small 
compared to the Fermi energy p^/2'm. In the absence of 
electron-electron interactions, functions G^ and G^ are 
independent of the non-equilibrium state of the system. 
In the Fourier representation, they are given by 
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Here £,p — {p^ — pji^)/2m is the kinetic energy counted 
from the equilibrium chemical potential, Ap — ap is the 
energy of the spin-orbit splitting, and fjp — i) ■ p/p is the 
projection of the spin operator rj onto the direction of the 
electron momentum. The Keldysh function G^ satisfies 
the equation 
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It is now customary to apply the Wigner transformation 
to Eq. (O, i.e. the Fourier transform with respect to the 
relative time and space arguments. 
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here t± = tztSt and x± = x ± Sx/2. In the semiclassical 
approximation, the Wigner transform of the right-hand 
side of Eq. (O can be replaced by a product of the Wigner 
transforms of S and G: 
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where V = V + eE d^ , and 



is the density matrix of electrons with the energy e. The 
total number of particles and their total spin can be ex- 
pressed via pe as follows. 



N = Tri^ depe, 



S = iTr. 
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In the limit r — > cx) the equation Q reduces to the bal- 
listic equation of Ref. [l^- Note, however, that the func- 
tion ^pe is not a distribution function in the conventional 
sense, since it depends on both energy and momentum. 
A stationary solution to the quantum kinetic equation 

is of the form ^p^ = ^^(Gpg — Gp^), where A^ is an 
arbitrary scalar function of the electron energy s. This 
solution represents the state in which the charge density 
is uniform, and spin density is zero. In a non-equilibrium 
state with the characteristic length scales of the spin and 
charge densities exceeding the electron mean free path 

1 = vpT, the distribution gp^ relaxes slowly to equilib- 
rium. To describe this relaxation, we derive the equation 
for the density matrix /5e(r,t). It is useful to move small 
gradient terms to the right hand of the kinetic equation 
(0 , so that its left hand side describes fast relaxation to 
the local equilibrium distribution: 



where 
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Small anisotropic deviations from local equilibrium are 
due to the gradient term ICpJ in the kinetic equation 
which can be treated perturbatively. The solution to 
Eq. H10|l can formally be written (in the Fourier repre- 
sentation with respect to time) as, 
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where 57 — oj + iIt. In a zeroth order, one can neglect the 
gradient term ICpe altogether, so that Eq. ifT^ gives the 
distribution gpj in terms of the density matrix p^. In the 
first order, we substitute the obtained expression for gp^ 
in the gradient term JCpe to obtain an improved expres- 
sion for the distribution function, gpe ■ This procedure is 
then to be repeated to the necessary order. 
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Integrating the second-order approximation over the mo- 
mentum p, one arrives at the diffusion equation for the 



density matrix p^. In a quasistationary regime [lot ^ 1) 
the equation takes the following form: 
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(14) 
The first two terms in this equation describe spin and 
charge diffusion with D = VpT/2 being the conventional 
diffusion constant, and vp — pp/f^ the Fermi velocity. 
The third term describes a spin precession due to the 
drift velocity, and the fourth term describes the cou- 
pling between charge and spin. The right hand side of 
Eq. (|14|l describes spin relaxation due to the Dyakonov- 
Perel mechanism [l^ . The coefficients of the Dyakonov- 
Perel spin-relaxation, spin-density coupling and spin- 
precession are, respectively, 
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where A = Ap^ , and the dimensionless parameter ( = 
Ar describes relative strength of spin-orbit coupling and 
disorder scattering. In deriving Eq. 115(1 . we assumed 
that the spin-orbit splitting is small compared to Fermi 
energy (A <C Ep), while the parameter C — At is arbi- 
trary. (Physically, ( represents the angle of spin preces- 
sion between two consecutive collisions.) In the case of 
weak spin-orbit coupling or a very clean sample (C <C 1), 
the Dyakonov-Perel relaxation time is large compared to 
the elastic mean free time Tg ^ t/C^ 3> t and the char- 
acteristic spin relaxation length ^fDrl is large compared 
to the mean free path. The spin dynamics is thus slow 
both in space and time and Eq. (|14|l has a meaning of a 
real diffusion equation for the coupled density and spin 
degrees of freedom. Note that in the limit C ^ 1 the 
spin-density coupling (i3-term) in the equation l|14|) dif- 
fers from the corresponding term in Ref. [ill ^^ ^ result of 
incorrect summation of a diffusion ladder in Ref. [Il| . We 
show below that the strength of this coupling is crucial 
for the magnitude of the spin-Hall effect. 

In the opposite limit, C ^ 1: spin-relaxation is fast, 
Ts ^ T, and occurs on a length scale of the mean free 
path /, i.e. locally as compared to the system size L ^ 
I. Spin relaxation dynamics (e.g. propagation of a spin- 
polarized injected beam) is therefore beyond the reach 
of the diffusion equation and must be studied with the 
kinetic equation iQ. However, Eq. H14|) can still be used 
to study a steady state in which spin polarization changes 
slowly on a scale of I (which will be the case for spin- 
Hall conductivity, see below) . One then has to retain the 
terms describing density diffusion, spin relaxation and 
spin-density coupling. In the vector basis. 
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the equations H14|) are reduced to. 
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Total density and spin polarization are expressed in this 
basis as: N = v J den,;, and S — v J des^. 

Spin- accumulation. We now apply the spin diffusion 
equation l(14() to analyze spin accumulation in a finite- 
size sample of the length L contacted by two ideal un- 
polarized metallic leads. The sample is infinite in the 
transverse direction so that Pe{x) depends on the longi- 
tudinal coordinate x only. Note that the electric field 
in the sample enters Eq. (14) only via V = V + eE9e 
and therefore can be eliminated by shifting the energy as 
e ^- e -I- eEx. Thus, the electric field may be treated via 
the boundary conditions. 
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where V — EL is the voltage bias between the two leads, 
and Fg is the equilibrium Fermi-Dirac distribution. Sub- 
stituting the expansion (|16|) into Eq. (|14|l we observe that 
s^ = si = 0. The other two equations yield. 
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Note that the B-term in the equation for n^ leads to small 
corrections, ~ a'^/vp, which must be neglected in the 
considered approximation. The solution to the second of 
Eqs. (dnj yields. 
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cosh [7 — x/Ls] 
cosh 7 
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where 7 = L/2Ls is the dimensionless spin-flip rate, and 
Ecs ~ V/L. For an infinite system, 7 — > cx3, the spin 
accumulation 12011 agrees with the previous calculation 
by Edelstein jl4|. 

Spin- current. The spin current, as defined in the in- 
troduction, is found from the Keldysh Green's function. 
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The function G^ can be expressed via the density p^ with 
the help of the equation, G^ = G^T,^G^, which follows 
from Dyson's equation Q. After simple transformations. 



fk = -o Tr<T^(V^-Vfc) 



xp,(y)G^(y-x') + -Qfc.A^. 
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Keeping now in the integrand only the zero and first- 
order terms in the expansion of p^ over y — x, we arrive at 
the final expression for the non-equilibrium spin current 
in terms of the density and spin distribution functions, 
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dependent, spin polarization is retarded with respect to 
the field, due to the finite spin relaxation time. As 
a result, the spin polarization contribution in Eq. II23I) 
does not exactly cancel the electric field contribution, 
and spin-Hall conductivity becomes non-zero. Solving 
Eq. H1U|I for the homogeneous infinite system and gen- 
eralizing Eq. (|22l) for a time-dependent state, we find. 



FIG. 1; (color online) In a spin-Hall bar setup, the electric 
current is driven through 2DES contacted by the metallic 
leads connected to a voltage source. The electric field E^ 
creates an in-plane spin polarization Sy in the bulk. Spin 
currents jy are running in the vicinity of the contacts while 
vanishing in the bulk. Out-of-plane polarization Sz is accu- 
mulated at the sample corners. 



Here Ecff = E — 'VN /2ev is the gradient of the electro- 
chemical potential including both the electric field and 
the gradient of electron density. 

Substituting Eq. ^^ into Eq. ((^SJ we observe that the 
two contributions to j], cancel each other in the bulk of a 
sample. Therefore, the dc spin current vanishes indepen- 
dent of the relative strength of disorder and spin-orbit 
interaction. Thus, we generalize the result by Inoue et 
al. [9j to arbitrary values of (. (The finite value of the 
spin-Hall conductivity obtained in Ref. [12| was due to 
mishandling of the electric field vertex in the calculation 
of the spin-Hall conductivity.) 

However, near the contacts where the spin polarization 
deviates from its bulk value, the spin current is non-zero. 
Using the expression (^0)) in Eq. if^ . we find that the 
spin current near the contacts decays as (C «C 1), 
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For a sample of finite width, this spin current should 
lead to a spin accumulation within a distance Ls of the 
corners of the sample, with a component of S along ±z, 
as illustrated in Fig. 1. 

Note, that for a non-uniform system in thermal equilib- 
rium, where Eeff — 0, the spin density given by Eq. l(T^ 
is zero, as well as the spin current. Small equilibrium spin 
currents [Sj, proportional to {a/vp)^, are beyond the ap- 
proximation used when deriving Eq. H23|) . Our derivation 
of the diffusion equation H14|) and the spin-ciurent H23|) 
relies on the approximation Q that neglects contribu- 
tions from diagrams with crossed impurity lines (ladder 
approximation). This is usually justified provided that 
EpT 3> 1. In an infinite system the result H23|l is equiv- 
alent to a calculation within the Kubo formalism with 
the first term representing a single-loop contribution and 
the second term originating from the ladder impurity di- 
agrams. 

To reconcile our result for spin current with the pre- 
dictions of Refs. 13 , it is helpful to consider the ac spin- 
Hall effect [9J. When the applied electric field is time- 



For low frequencies, loTs < 1, the spin-Hall conductivity 
remains small, asH ~ —itor. When the frequency ex- 
ceeds the spin relaxation rate {ojTs > 1), asH reaches its 
maximum value er / (AiTTg) . For clean samples, this is the 
universal value c/Stt predicted in Refs. |2|, while for dirty 
samples (C ^ 1) the maximum value of the spin-Hall con- 
ductivity remains strongly suppressed, asH ~ eC^/(27r). 
To conclude, we derived quantum kinetic equation for 
2D electrons in the presence of spin-orbit coupling and 
short-range potential scattering. We proved that the dc 
spin Hall effect disappears in a bulk sample, and we com- 
puted the spin accumulation in a finite-size system for a 
wide range of parameters. This work was supported by 
NSF grants PHY-01-17795 and DMR02-33773. 
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